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Abstract. We study, by means of mean field calculations and Monte Carlo simulations of a lattice-gas 
model, the distribution of adhesion sites of a bilayer membrane and a supporting flat surface. Our model 
accounts for the many-body character of the attractive interactions between adhesion points induced by 
the membrane thermal fluctuations. We show that while the fluctuation-mediated interactions alone are 
not sufficient to allow the formation of aggregation domains, they greatly reduce the strength of the direct 
interactions required to facilitate cluster formation. Specifically, for adhesion molecules interacting via a 
short range attractive potential, the strength of the direct interactions required for aggregation is reduced 
by about a factor of two to below the thermal energy fcsT. 
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1 Introduction 

Adhesion between two membranes or between a mem- 
brane and another surface is an important topic for its 
ubiquitous occurrence in biological and biophysical pro- 
cesses. This process, during which two interfaces attract 
each other, can in principle be facilitated by non-specific 
attractive interactions (e.g.. Coulomb and van der Waals 
interactions). Cell adhesion, however, is usually caused by 
highly specific receptor molecules located at the outside 
of the plasma membrane of the cell, that can bind to 
specific ligands on the opposite surface [HH]. Typically, 
the area density of receptors is rather low which does not 
lead to efficient bio-adhesion. However, when facing a sur- 
face with enough ligands, the receptors may cluster into 
highly concentrated adhesion domains to establish much 
stronger binding [S'lT . Formation of adhesion clusters oc- 
curs in many biological processes [5^ , including the binding 
of white blood cells to pathogens |5], cadherin- mediated 
adhesion of neighboring cells [7] , and focal adhesion of cells 
to the extracellular matrix [5]. Much insight into these 
bio-adhesion processes has been gained from experimental 
studies of biomimetic membranes with receptors molecules 
that interact with surfaces covered by ligands [Sl ITUllTTlfT^ 

mm- 

Generally speaking, adhesion induced domain forma- 
tion requires some attractive intermolecular interactions 
between the receptor-ligand pairs. These interactions in- 
clude both direct and membrane-mediated contributions. 
The former are typically described by pairwise interac- 
tions which are infinitely repulsive at very small molecular 
separations and attractive at somewhat longer (but still 
finite) distances [TS] . The effect of the direct interactions 



between adhesion bonds can be studied in the framework 
of the thoroughly researched lattice gas model JJ, . Specif- 
ically, there exists a critical value Cc of the strength of the 
attractive part of the intermolecular pair potential above 
which the system may phase separate into domains with 
high and low concentrations of adhesion bonds. 

Much less is known about the indirect interactions be- 
tween adhesion sites which are mediated by the mem- 
brane thermal fluctuations. The entropic origin of these 
interactions can be easily understood as follows: Consider 
two adhesion bonds between two membranes or between 
a membrane and a surface. The adhesion points restrict 
the thermal height fluctuations of the membrane in their 
vicinity. This entropy loss can be minimized if the two ad- 
hesion bonds are brought to the same place, in which case 
the membrane becomes pinned at only one point and not 
two. The membrane fluctuations, thus, induce an attrac- 
tive potential of mean force between the adhesion points. 

In a previous publication, we analyzed the membrane 
mediated interactions between two adhesion bonds of a bi- 
layer membrane and a supporting surface [T7]. We found 
that for "point-like" adhesion molecules whose size is com- 
parable or smaller than the thickness of the membrane 
(Z ~ 4 — 5 nm), the potential of mean force is an infinitely 
long range attractive potential that grows logarithmically 
with the pair distance r: 



C/(r) = 2fcBrin(y) 



(1) 



Eq.([T]) holds for tensionless membranes, which will be the 
case discussed in the following. When surface tension is ap- 
plied, this form holds for pair separations r ^ ^o- = -^/k/ct, 
where k and a denote the membrane bending rigidity and 
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surface tension, respectively. For r :^ £_cr,U{r) is decreased 
by a factor of 2 ^7\. We leave the discussion in stressed 
membranes to a future publication. 

In this paper we consider the same system as in ref . |17| , 
but instead of two adhesion points we look at a membrane 
which is pinned to a flat impenetrable surface at multiple 
sites. To analyze the aggregation behavior of the adhe- 
sion sites we first need to generalize eq. ([T]) and write down 
the fiuctuation-induced interaction energy as a function of 
the coordinates of the adhesion sites C/(ri, r2, ra, . . . , r^)- 
This is a non-trivial problem since the fluctuation-mediated 
interaction is a many-body potential which cannot be ex- 
pressed as the sum of two body terms of the form in 
Eq.([T]). The many-body nature of U{ri,r2, r^,. . . , r]\r) is 
best illustrated by an example: Consider a cluster of two 
adhesion points located close to each other at ri ~ r2, 
and a third distant adhesion point located at r^. Hav- 
ing a single adhesion point at ri or r2 instead of the 
two-point cluster will not result any change in the spec- 
trum of membrane thermal fluctuations. Therefore, the 
third point is attracted to the two-point cluster by the 
same potential of mean force ([1]) to which it is attracted 
to one adhesion point, and not by a potential which is 
twice larger than potential ([T|), which would be the case 
if C/(ri,r2,r3, . . . ,rjv) is the sum of pair interactions. 
What is the exact form of C/(ri, r2, ra, . . . , r]\[) is still an 
open question that needs to be resolved for the fluctu- 
ation induced domain formation to be understood. Sev- 
eral approximations to this problem which avoid direct 
many-body calculations have been proposed. Weikl and 
Lipowsky introduced a mean field theory in which the ef- 
fect of the pinning points is represented by a homogeneous 
attractive interaction between the fluctuating membrane 
and the underlying surface (or between the membrane and 
another membrane) (TB]. They concluded that the homo- 
geneous fluctuation induced potential alone cannot facil- 
itate the formation of adhesion zones, but it greatly re- 
duces Cc, the critical strength of the direct interactions 
between the adhesion bonds above which the formation of 
adhesion clusters is possible. A very similar conclusion has 
been recently reached by Speck et al. using rigorous sta- 
tistical mechanical methods 19 . However, in their model 
the hard wall interaction between the membrane and the 
surface has been replaced by a harmonic confining poten- 
tial. 

In this work we present a more accurate approach to 
the problem which employs a non-additive many body po- 
tential U{ri,r2, ra, . . . , rjv). The derivation of the poten- 
tial [/(t"!, r2, ra, . . . , r^), which is based on our previous 
statistical mechanical studies of the membrane thermal 
fluctuations with one [5(T and two [T7] adhesion points, 
is presented in section [51 The idea is somewhat different 
than the one previously used by others. As in refs. P^ITO] . 
we integrate out the membrane degrees of freedom and 
map the problem onto the lattice gas model where the oc- 
cupied sites represent the adhesion bonds. But unlike in 
refs. [T81IT9] . the interaction energy f/(ri, r2, ra, . . . , rjy) 
is not expressed as the sum of two body interactions be- 
tween nearest neighbor lattice sites. Instead, it is calcu- 



Fig. 1. Schematic of the system under investigation, consisting 
of a membrane which is pinned to a flat impenetrable surface 
at multiple adhesion points. The adhesion points can diffuse 
freely on the surface, and they may cluster to increase the 
conformational entropy of the fluctuating membrane. 



lated from the empty sites that represent the unpinned 
sections of the membrane. The energy assigned with each 
empty site represents the loss of entropy due to the re- 
duction of the membrane fluctuations in the area around 
that site. We argue that the extend by which the mem- 
brane fluctuations are locally restricted depends mainly 
of the distance to the closest pinning point. Therefore, 
the calculation of U{ri,r2, r^, . . . , r^) for a given distri- 
bution of pinning sites involves the division of the lattice 
into Voronoi cells and summing their contributions to the 
free energy. In section [3] we present a mean field analy- 
sis of our model, and in section S] we present the results 
of Monte Carlo lattice simulations. We find, in agreement 
with other theoretical studies, that while the fiuctuation- 
mediated interactions alone are not sufficient to allow the 
formation of aggregation clusters, they greatly reduce the 
strength of the residual (direct) interactions between ad- 
hesion points which is required to facilitate cluster forma- 
tion. More specifically, we find that the strength of the 
short range cohesive energy that allows condensation is 
reduced by about a factor of two and falls below the ther- 
mal energy fc^T. In section [5] we summarize and discuss 
our results. 



2 The many body fluctuation-mediated 
potential 

We consider the system shown schematically in fig. [1] con- 
sisting of a fiuctuating membrane of linear size L which is 
pinned to a fiat impenetrable surface at several sites. The 
free energy cost of a single adhesion point is given by pTl 



Fi 



.Tlni^ 



2kBT\n{ J 



(2) 



This result has been derived in ref. (20j by noting that 
pinning the membrane to the surface at one point does not 
modify the membrane spectrum of thermal fiuctuations. 
It does, however, eliminates the membrane translational 
degree of freedom by enforcing the global minimum of 
the membrane height function to be located at the point 
of contact with the surface. A different interpretation to 
eq. ^ is that between the membrane and the surface there 
is an effective interaction free energy per unit area of the 
form 

nr) ^ ^, (3) 

that represents the entropy cost due to the reduction of 
the membrane fluctuations at distance r from the adhesion 
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site [T7]. The attachment free energy is then derived by 
integrating the free energy density U{r) over the projected 
area of the membrane 



Fi 



2TrrV{r)dr. 



(4) 



The free energy density ([3]) can hkewise be employed 
to derive the fluctuation mediated pair potential ([T]) . How- 
ever, when we have two adhesion points, the distance r in 
eq.® should be the distance of the unit membrane area 
to the closest adhesion point. With this choice of r, let us 
now consider a square membrane {—L/2 < x,y < L/2) 
with two adhesion points located at {x,y) — (±ro/2,0). 
The attachment free energy is calculated by integrating 
V{r) over the projected area of the membrane (excluding 
a region of size / around the adhesion points): 



F2 = 4 



L/2 



dy 



('■0-0/2 



dx ■ 



knT 



TT y^ + {x - ro/2)" 



L/2 



dx ■ 



kuT 



{ro+l)/2 TT y^ + {x - ro/2y 



(5) 



Integrating over y yields, 



F, 



^ksT 



L/2 



('■0-0/2 



dx 



\x-ro/2\ 



■ tan 



2\x~ro/2\ 



(ro-0/2 



dx 

^ro/2 



tan 



L 



2\x-ro/2\ 



(6) 



Assuming that I < vq <S^ L, the inverse tangent function 
in eq.® can be approximated by the constant value of 
7r/2 over most of the integration range (except near the 
boundaries of the membrane |a::| ^ L/2 which, neverthe- 
less, does not influence the dependence of the result on tq). 
With this approximation, one gets 



F2iro,L) ~ 2A:BTln - + 2fcsTln ^ 



L 



ro 



FiiL) + U{ro). 



(7) 



The first term in eq.© is the free energy cost of a sin- 
gle adhesion site ^, which is the expected value when 
the two adhesion points coincide (rg ~ ^) to form a single 
cluster. The second term, which represents the additional 
free energy cost associated with the separation of the ad- 
hesion points, is identified as the fiuctuation induced pair 
potential, in agreement with eq.(IT]). 

We now wish to generalize eqs.© and ([7]), and pro- 
pose that the many-body fiuctuation-mediated potential 
is given by 

C/(ri,r2,r3,...,rjv) = Fn {ri,r2,r3, . . . ,rN,L)-Fi{L). 

(8) 
The attachment free energy Fn of N adhesion points is 
calculated by integrating the free energy density 



F 



N 



kuT 



7rd2(r,ri,r2,r3,...,rAr) 



dr, 



(9) 



where d is the distance of the membrane unit area to the 
nearest adhesion point 



d = min 

i=l...N 



\r-ri\)^ 



(10) 



and the integration is carried over the projected area of 
the membrane except for small regions of size I near each 
adhesion point. What we essentially argue here is that 
the extent by which the membrane thermal fluctuations 
are limited at each spot depends almost exclusively on 
the distance to the nearest adhesion point. All the other 
adhesion points are effectively screened. This suggestion 
is supported not only by the above example, eqs.®-©, 
but also by our observation from recent computer simula- 
tions presented in ref. [17 . In that paper, we determined 
the pair potential (HJ by using Monte Carlo simulations of 
a coarse-grained bilayer model with two adhesion points. 
Our simulations results were in excellent agreement with 
eq. (dl despite of the fact that each adhesion point inter- 
acts not only with the other adhesion point but also with 
its infinite array of periodic images. It is unlikely that 
the periodic images have such a negligible contribution on 
U(r) for all values of r unless they are simply screened. 
In other words, the local fluctuation behavior of the mem- 
brane is governed by the distance to the nearest adhesion 
point while the spatial distribution of the other, more dis- 
tant, points is irrelevant. Thus, calculating Fjsi for a given 
distribution of adhesion points involves the construction 
of the 2D Voronoi diagram of the conflguration, and sum- 
mation of the free energy contributions coming from each 
Voronoi cell, Si: 



N 



F, 



N 



E ^^^ 



kuT 



(11) 



where r is the distance to the adhesion point which is 
located inside the Voronoi cell. 



3 Mean field theory 

To study the formation of adhesion clusters in supported 
membranes, we consider a lattice of Ns sites (with lattice 
spacing set equal to I (the cut-off microscopic length scale) 
of which N < Ns sites are occupied by adhesion points. 
For each spatial configuration of the system, the energy is 
written as the sum of two terms 



E = El 



LG 



"Af- 



(12) 



For the first term in ea. p2l) representing the direct inter- 
action between the adhesion points, we take the lattice 
gas energy with nearest neighbor interactions 



(id) 



:/' 



(13) 



where Si = 1 for an occupied site, s^ = for a vacant site, 
and e is the site-site interaction energy. The second term 
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in eq. (J12p , which represents the attachment free energy, is 
given by the discrete form of eq. ([TT|) 






(14) 



where r is the distance from a vacant lattice site to the 
nearest occupied lattice site. 

The phase behavior of the system can be analyzed 
through mean field theory. Let us assume that the ad- 
hesion points form Nc < N clusters. The free energy of 
system includes three contributions: (i) the mixing entropy 
of the adhesion clusters, (ii) the energy Ei^q of the direct 
interactions between the adhesion points, and (iii) the at- 
tachment free energy F^. The first free energy contribu- 
tion is given by 



F^ 



ksT 



N, 






1 



^=(fi. <-• 



where c is the second virial coefficient. On average, each 
cluster consists of (N/Nc) adhesion points; and if we as- 
sume that it has a roughly circular shape than c ~ 4(iV/7Vc). 
Denoting the number densities of the adhesion points by 
cj) = N/Ns, and of the clusters by (j)* = N^/Ns < 4>, the 
free energy of mixing per lattice site is given by 



^n- 



NskeT 



=.0* [In (</)*) -l]-t-2(/)(^* 



(16) 



The second contribution to the free energy is due to 
the direct interactions between the adhesion points. The 
ground state of the interaction energy ii^LG is achieved 
when a single circular adhesion domain with minimal sur- 
face is formed. If we set the ground state as the reference 
energy, the energy of an ensemble of clusters can be es- 
timated as being proportional to the total length of the 
domain boundaries. For Nc circular clusters of size [N/Nc) 
we have 



^LG ^,Nc N_ 
NskeT Ns V Nc 



(17) 



where A, the associated dimensionless line tension, is pro- 
portional to the interaction energy e in ea. (fT3|) and S, the 
mean number of nearest-neighbor vacant sites per occu- 
pied site on the boundary of a cluster 



X = Be. 



(18) 



The sum of free energy contributions P^ and p7)) con- 
stitutes the total free energy density (per lattice site) of a 
2D lattice gas of clusters: 



Flc 



NsksT 



^ (j)* \n{(j)*) - (j)* + 2(j)(j)* + Xy 



(19) 



The third contribution of the attachment free energy 
can be estimated as follows. The clusters form Nc Voronoi 
cells, each of which has on average an area of Avor = 
{Ns/Nc)P ■ The attachment free energy of each Voroni cell 



is given an equation similar to eq.([2]) for the attachment 
free energy of one adhesion point, but with v4vor instead 
of the total membrane area L^ . Thus 



F, 



N 



Nc 



fcsTln 



and the attachment free energy density is given by 
Fn 



NskeT 



= -rin(0*). 



(20) 



(21) 



which eliminates the first term in the lattice-gas free en- 
ergy density [ea. lfTO]) ]. yielding 



F 



F 



LG 



F, 



N 



NsksT NskeT NskeT 

(22) 
We consider a low density of adhesion sites 0^1, 
which also implies a low number density of adhesion clus- 
ters since (jf < (f). By minimizing the free energy density 
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Fig. 2. The phase diagram of the adhesion sites calculated 
within the mean field approximation. (A) Eq. (|23P for the 
lattice-gas model. (B) Eq.((24]) for adhesion points of fluctuat- 
ing membranes. Ai and A2 represent the first-order transition 
and spinodal lines, respectively. 
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we obtain the equilibrium value of the (jf for the lattice- 
gas problem [ea. (fT9)) ] and for the adhesion points of a fluc- 
tuating supported membrane [eg.ip^]. In both cases, the 
system undergoes a first order phase transition at Ai (0) 
from the gas phase (0* = </>) to a condensed phase consist- 
ing of only a few clusters (0* ^ 0). Also, in both cases the 
free energy reaches a maximum at intermediate densities 
(0 < 0* < 0). This free energy barrier for condensation 
disappears at the spinodal point \2{4>) > Ai (</>). For the 
lattice-gas problem we find 



Ai'^ = 1 - 20 - ln((/)) 



A^^ = -40-21n(<^), 



(23) 



while for the adhesion points of fiuctuating membranes we 
have 



Ai = 1 - 20 

Aa = 2-40 = 2Ai. 



(24) 



The results of eas. (P5)) and ([M]) are summarized in fig.[5]4. 
and B, respectively. The important points in the results 
are that: (i) Ai > 0, which means that the fiuctuation 
induced interactions alone are not sufficient to induce ag- 
gregation of adhesion domains, but (ii) they greatly reduce 
the strength of the direct interactions required to facilitate 
cluster formation since Ai < A^^ (and also A2 < Aa^). 
In the following section we support these conclusions with 
MC simulations. We show that for adhesion points of fiuc- 
tuating membranes, the site-site cohesive energy e for the 
onset of aggregation falls below the thermal energy fc^T. 



4 Computer simulations 

To further investigate the aggregation behavior of adhe- 
sion points, we performed Monte Carlo (MC) simulations 
of our lattice-gas model with the total configurational en- 
ergy given by the sum of direct [ea. (ll3l) ] and fluctuation- 
induced [ea.([Hl)] interactions. We used a 120 x 138 trian- 
gular lattice (that has an aspect ratio very close to 1) with 
periodic boundary conditions. We simulated the system at 
two different densities (f> — N/Ns = 0.05 and = 0.1, and 
for various values of e ranging from to 3 fcsT. For com- 
parison, we also simulated the standard lattice-gas model 
[for which the configurational energy is given by eq. (jl3p . 
without the fluctuation-mediated free energy ea. (fT4|) ]. For 
each density and for each value of e, we performed 8-16 
independent runs starting from different initial configura- 
tions where the points are either randomly distributed on 
the lattice (as in fig. [3J\) or put in a single cluster (see 
fig. [33) • The system was then equilibrated until the dis- 
tribution of points in all the independent runs look similar 
(see e.g., figEP us. |3P, andOil ws. [3F). Equilibrium time 
for the different samples ranges from 3.6 x 10^ to 10^ MC 
time units, where each MC time unit consists of N single 
particle move attempts. For the adhesion points problem, 
each particle was displaced to a randomly chosen nearest 
neighbor lattice site, which enabled us to employ an effi- 
cient algorithm to update the Voronoi diagram needed for 
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Fig. 3. Initial configurations of the simulations in which (A) 
the sites are randomly distributed on the lattice, and (B) put 
in a single compact cluster. Representative equiUbrium con- 
figurations of (C-D) our model [Eqs.([l3l) and ^] and (E-F) 
the standard lattice-gas model [Ea. H13p only] for <j) = 0.1 and 
e = IksT. Configurations (C) and (E) evolved from the initial 
state (A), while (D) and (F) evolved from (B). 



calculating the fluctuation-mediated free energy (J14p . For 
the standard lattice-gas model, each move attempt con- 
sisted of randomly selecting a particle and moving it to 
the nearest vacant point in a randomly chosen lattice di- 
rection [32] . After the first stage of equilibration, the sim- 
ulations were continued for 3 x 10"" MC time units during 
which data was collected every third MC time unit. 

To examine the occurrence of a phase transition from a 
gas to a condensed phase, we measured the average num- 
ber of clusters in the system (where a cluster is defined as 
a set of neighboring occupied sites), and the mean value 
of the energy of direct interactions between sites, (-Elg) 
[see ea. (|13p ]. Our results are summarized in fig. 2)4. (for 
(p = 0.05) and B (for </> = 0.1). For each </> we measured 
these quantities both for the standard lattice-gas model 
(open symbols and dash-dotted lines in fig. S]) and for 
adhesion points which also interact via the fiuctuation- 
mediated free energy ea. p^ (solid symbols and solid lines 
in fig. 21). The number of clusters is denoted by squares 
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3000 



800 




Fig. 4. Left y-axis: The energy of direct interactions between 
sites, (-Elg) [Ea. H13|l ]. as a function of e, for = 0.05 (A) and 
(j) = 0.1 (B). Sohd circles - results for our model for adhesion 
points. Open circles - results for the standard lattice-gas model. 
Right j/-axis: The number of clusters as a function of e, for 
4> = 0.05 (A) and (fi — 0.1 (B). Solid squares - results for 
our model for adhesion points. Open squares - results for the 
standard lattice-gas model. 



(values on the right y-axis of the figures), while (£^lg) is 
represented by circles (values on the left y-axis). 

The gas phase is characterized by a large number of 
small clusters, some of which may be of the size of a single 
site. Furthermore, since each occupied site has a relatively 
small number of neighboring occupied sites, the mean con- 
figurational energy (—-Elg) is relatively low. Conversely, 
when the sites form large clusters in the condensed phase, 
{—Ei,g) is high, and the total number of clusters decreases 
(and in many cases, especially for large values of e, we sim- 
ply observe only a single cluster in our system). Fig. 2] ex- 



hibits an abrupt, clearly first-order, transition from a gas 
phase with a large number of clusters and small (— Slg) 
to a condensed state with a small number of clusters and 
large (—-Elg)- The estimated values of e at the transition 
are (see vertical lines in fig.S]): et ~ O.T/cbT {4> = 0.05) and 
et ~ O.GSfcsT (0 = 0.1). In comparison (see also fig.|ll), for 
the standard lattice-gas model, the transition values are 



roughly twice larger than these values: e 



LG 



l.AbksT 



{(f) = 0.05) and e^G ^ LSkeT (0 = 0.1). Fig. [3p-F 
exhibits typical equilibrium configurations of the system 
at = 0.1 and e = Ifc^T. For the lattice-gas model 
ej"*^ > IksT, and at equilibrium the system is in the gas 
phase (figs. 135 and F). When the fluctuation-induced in- 
teractions eg. pi)) are introduced, Cf falls below IfcsT and 
the system is in the condensed phase where most of the 
particles belong to one large cluster (figs.[3p and D). 

Our computational results which show that the fluc- 
tuation mediated interactions reduce the strength of et, 
are in a qualitative agreement with the mean field theory 
prediction (section |3]). To make a quantitative compari- 
son between the theory and the simulations, one needs to 
estimate the parameter B appearing in ea. p^ . Several 
reasons make such an estimation difficult and inaccurate: 
First, our non-standard mean field theory is based on the 
assumption that the clusters are circular and roughly have 
the same size, which is quite a crude approximation. Sec- 
ond, tracing the precise location of Ct in fig. |4]) is largely 
inaccurate because of the finite size of the system that 
makes the transition look like a crossover. To reduce the 
large uncertainties associated with the determination of 
et, one can look at the difference between the value of this 
quantity in our model and for the standard lattice-gas 
model. Using 



A;^^ - Ai = B 






(25) 



for (/) ~ 0.1, we find B ^ 3.5 which for the triangular 
lattice with 6 nearest-neighbors means that adhesion point 
residing on the edge of a cluster loose slightly more than 
half of their neighbors. 



5 Summary and discussion 

In this paper we studied the aggregation behavior of ad- 
hesion points between a fiuctuating membrane and a sup- 
porting surface. We demonstrated that the problem can be 
mapped onto a lattice-gas model with two types of molec- 
ular interactions: (i) direct site-site pair interactions and 
(ii) Casimir-like interactions which are mediated by the 
membrane thermal fluctuations. The fluctuation-mediated 
interactions, which are inherently of many-body character, 
are calculated in our model by summing over the vacant 
rather than the occupied sites of the lattice. Each vacant 
site represents a small unit area of the fluctuating mem- 
brane, and the fluctuation-mediated potential expresses 
the local free energy cost due to the restriction imposed 
by the adhesion points on the membrane thermal fluctu- 
ations. This free energy cost depends mainly on the dis- 
tance between the vacant sites and the nearest occupied 
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site. Therefore, such a many-body potential is calculated min(l,exp(— zlS/fcsT)), where AE is the energy change 

by determining the Voronoi diagram for each lattice con- caused by the move attempt, 

figuration, which can be quite easily implemented in MC 
simulations. 

We used mean field calculations and MC computer 
simulations to investigate the phase behavior of a lattice- 
gas of adhesion sites at low densities. We showed that 
upon increasing the strength of the site-site interactions 
e, the system undergoes a first-order phase transition into 
a condensed state. The fluctuation-induced interactions 
lower the value of e at the phase transition to below the 
thermal energy fc^T. This result suggests that fluctuation- 
mediated effects play a central role in the formation of ad- 
hesion domains in biomimetic and biological membranes. 

This work was supported by the Israel Science foun- 
dation (Grant Number 946/08). 
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